AN EFFICIENT METHOD OF FINDING A PUISEUX 
EXPANSION OF A PARAMETRIC SINGULARITY 
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Abstract. We present a new and very efficient method of determining 
a topological type of a cuspidal plane curve singularity parametrised 
locally by (x(t),y(t)). The approach relies on defining inductively a 
sequence of functions Pk which are polynomials in x, y and derivatives 
of x and y. 



Let us be given a germ of a parametric curve in C 2 given by 
x(t) = a t p + ait p+1 + . . . , 

(1) y(t) = M 9 + M 9+1 + ..-. 

Here anfrn 7^ 0, q > p > 1 and both series on the right hand side are assumed 
to be convergent in some neighbourhood of t = 0, In order to determine 
the topological type of the singularity JT]) we need (see [ BKl IEN] ) to find its 
characteristic sequence (p; qo, q%, . . . , q n ). To do this we write the Puiseux 
expansion of JT]) in the form 

(2) y = c q x q / p + c q+1 x( q+1 V p + .... 

Then we define qo = min{s > 0: c s / A p \ s}, and inductively po = p, 
p k = gcd(pjfc_i, gjfe-i). If Pk ¥= 1 we define q k = min{s > 0: c s / A p k \ s}; 
if Pk = 1 we stop the procedure. As p k > Pk+i this procedure will definitely 
stop after a finite number of steps. So let n be such that p n +i = 1, then 
(p; qo, . . . , q n ) is the characteristic sequence. 

Therefore, in order to know the topological type of singularity we need to 
know the coefficients Ck, or, at any rate, to know which one are zero and which 
are not. Now passing from (JTJ) to (J2j) is, in general, a very involved procedure, 
which requires finding an explicit formula for t = d\x x l p + d<ix 2 l p + .... 
Moreover, the Puiseux coefficients c k are not very well defined, namely they 
are functions of some fractional powers of an. In particular if we study a 
deformation, in other words if we allow the coefficients of (pQ) to vary with 
some parameter s, and an(s) — > as s — > so, then c&(s) can easily get out 
of hand. Some of them can diverge to infinity, some of them can converge 
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but lose any topological meaning. We want to propose a remedy to both 
problems: we shall define some polynomials in y, x and their derivatives over 
t, which carry the same data about the singularity as Puiseux expansion and 
such that one can recover the Puiseux expansion from them in a very easy 
way. 

As a matter of fact, the definition of polynomials P k appears first in [Bo], 
but only after writing this article the author realised that the polynomi- 
als P k may be of wider use and allow fast and efficient computation of a 
characteristic sequence. 

We begin with putting 

P (t) = y and r = ord i=0 P (t) = q 

Like in [BZll Proof of Lemma 3.2], let us divide J2)) by x q l p and differen- 
tiate both sides with respect to t (it is convenient to write here x = ^| and 
so on). We get 

j^— = -c q+1 xx 1/p 1 + -c g+2 xx 2/p 1 + .... 

X1/P+ 1 p q p q 

So if 

(3) Pi = yx- -xy 



we will have 



P 



Pi = -Cq+lX^^X + -C q+2 X {q+2)/p X 



P p 

Let 

(4) ri = ardt= Pi-(p-l). 
Then we have an obvious lemma. 

Lemma 1. The coefficients c ro +i, . . . , c ri _i all vanish. Moreover the first 
not vanishing coefficient of Pi, i.e. the one at fi+(p-i) ; is equa i to ( n _ 

\ n/v+l 
r )c ri ■ a Q 

Now we define inductively the polynomials Pk and the numbers r^. The 
induction procedure is as follows. Assume that we have 

(5) Pk = • • . . . • razn^^^P + . . . 

p p p 

with 

(6) r fc = ord t=0 P fc -(2fe-l)(p-l). 

In ([5]) dots denote terms like (s — rk-i) . . . (s — ro)p~ k c s x 2k ~ 1 x s ^ p , which are 
too long to be written in a nice looking formula. 

Now we divide both sides of J5]) by x 2k ~ 1 x Tk ^ p , differentiate both sides 
with respect to t and then multiply them again by x 2k x T ' k ^ p+1 . We get 
r k+ \ - r k r k+ i - r k -i r k+ i - r . 2fc , , , , 
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where 

(7) r k+1 = ord t=0 P k+1 - (2k + l)(p - 1) 
and 

(8) P fc+ i = xxj t P k - {^x 2 + (2k - l)xxj P k . 

Similarly we have a lemma the proof of which is an easy computation. 

Lemma 2. The Puiseux coefficients c Tk _ x +i, . . . , c rfc _i all vanish. The co- 
efficient at t r fc+( 2fc_1 )(P _1 ) of P k , i.e. the first non-vanishing one, is equal 
to 

(r k - r fc _i)(r fc - r fc _ 2 ) ■ ■ ■ (r k - r )p k ~ 1 c rk a k/p+ ^ k ' 1 \ 

Thus, using the above method we can easily find the Puiseux expansion. 
The computations become simpler if there are many Puiseux terms that 
vanish. The algorithm, in a form easy to use, can be presented as follows. 

Algorithm. Given a parametric presentation as in ([TJ of a cuspidal singu- 
larity with aobo 0. 

Step 1. Set po = p and ro = q. Let p\ = gcd(poj<7o)- If Pi = 1 then the 

singularity is quasi-homogeneous and we stop. 
Step 2. Compute Pi as in ([3]) and T\ from (j3J|. Let p2 = gcd(pi,n). 
Step 3. If P2 = 1 then we stop. Otherwise put k = 1. 
Step 4. Compute P k +i and r k+ i from J8)) and j7| respectively. 
Step 5. Let p^+i = gcd(pfc,rfc) if Pk+i = 1 we stop, otherwise we increase k 

and go to Step 4. 

As an output we get a sequence of positive integers (f*o,ri,r2, . . . , r n ), such 
that the Puiseux expansion is 

y = c ro x r °/ p + c ri x ri//p + • • • + c rk x Tk ^ p + inessential Puiseux terms. 

Here c ro , c ri , . . . , c rk are all non-zero. The procedure will eventually stop 
unless the parametrisation ((T|) is not one to one. 

Of course, we can recover not only the characteristic sequence, but also 
the exact values of non- vanishing Puiseux coefficients. It is enought to look 
at the first non-vanishing coefficients of functions P k +\ and use second part 
of Lemma El 

Example. Consider a singularity parametrised by 
X (t) = t 12 + t 13 + |t 14 

y(t) = J* + -t 19 + -t 20 + -t 21 + —t 22 - —t 2 \ 

yy ' 2 14 14 1568 3136 

We write 

y = x Wl2 + Cl9X 19/12 + _ 
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But it is easy to compute that 

= 2025 35 _ 24975 36 
1 10976 43904 ' 
So r\ = 35 — 11 = 24. Thus Cig = C20 = C21 = C22 = C23 = and C24 / 0. 
Then 

= 2500875 59 
2 76832 

so T2 = 59 — 3 • 11 = 26, so C25 = and C26 7^ 0. Then a simple com- 
putation shows that P3 ~ t 82 , so C27 7^ and the Puiseux expansion is 
y = 2-3/2 _|_ CmX 2 _j_ c 26 x 13 / 6 + C27X 27 / 12 + . . ., so the characteristic sequence 
is (12; 18, 26, 27) (C24 is inessential). 

We can use the polynomials Pk in studying deformations. In fact, they 
behave very well if x(t) and y(t) vary smoothly. This approach is pursued 
in |Boj with quite a success in some partial cases. 
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